NASA TN D-909 


NASA TN D-909 



/// 


-ctv^ 


TECHNICAL NOTE 


D-909 


ANALYTICAL INVESTIGATION OF 

AN ADAPTIVE FLIGHT -CONTROL SYSTEM USING 

A SINUSOIDAL TEST SIGNAL 

By Jack E. Harris 

Langley Research Center 
Langley Field, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON June 1961 







IK 


L 

1 

% 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

TECHNICAL NOTE D-909 

ANALYTICAL INVESTIGATION OF 
AN ADAPTIVE FLIGHT- CONTROL SYSTEM USING 
A SINUSOIDAL TEST SIGNAL 
By Jack E. Harris 

SUMMARY 


An analytical study was made of an adaptive flight- control system 
which measures vehicle response to small- amplitude control-surface 
deflections produced by a sinusoidal test signal. Changes in the response 
to this signal are related to environmental changes, and the system is 
continuously altered to maintain this response equal to a preselected 
value. The system is suitable for use in high-performance aircraft and 
missiles and requires only the addition of a signal generator and a logic 
circuit consisting of a filter-rectifier network and a comparator- 
integrator network to a basic command- control system. Thus, it presents 
a relatively simple approach to the problem. 

The effects on system performance of variation in flight condition, 
system-gain level, test-signal frequency, and sensor location are included 
in the analysis. Longitudinal control of a high-performance research 
aircraft over flight conditions ranging from landing approach to a Mach 
number of 5*8 at an altitude of 150,000 feet, and longitudinal control 
of a four-stage solid-fuel missile including the first bending mode over 
the atmospheric portion of a launch trajectory constituted the basis for 
the analytical study. 

Results of an analog- computer study using time-varying coefficients 
are presented to compare the control obtained with the adaptive system 
with that obtained with a fixed-gain system during the atmospheric por- 
tion of a missile launch trajectory. 

The system has demonstrated an ability to maintain satisfactory 
vehicle control- system stability over wide ranges of environmental 
change. 
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INTRODUCTION 


The control problem associated with present-day aircraft and mis- 
siles has been greatly increased by the operational requirements of 
these vehicles. The wide ranges of Mach number and altitude, coupled 
with wide changes in center-of -gravity position (sometimes nearing 
20 percent of the vehicle length during the atmospheric portion of a 
missile trajectory), cause the stability and control characteristics to 
vary greatly. Often these characteristics vary so much that they are 
unsatisfactory during portions of the desired flight envelope. 

Automatic-control systems and stability-augmentation devices have 
been installed in most high-performance vehicles to provide satisfactory 
stability and control characteristics. Automatic longitudinal flight- g 

control systems in present use require variations of control-system 
parameters. These parameters are varied as programed functions of the 
flight conditions, such as airspeed, altitude, or time. For adequate 
programs, knowledge of the stability and control characteristics through- 
out the flight envelope is required. This knowledge in turn requires 
wind-tunnel tests, backed up by lengthy flight-test programs. 

A system capable of measuring its response, comparing this measured 
response with a desired standard, and modifying its parameters in a 
closed-loop fashion to obtain the desired response might greatly reduce 
the required flight-test prpgram and produce a more closely integrated 
control system. Such a system is called an adaptive control system. 

The potentialities of adaptive controls in handling the increased 
problems of control of aircraft have produced a great deal of interest. 

Most adaptive concepts can be summarized as follows: 

(1) Those that use the system response to a normal input, such as 
a pilot Input or external disturbance, to modify the system parameters. 

(2) Those that use the system response to an internal continuous 
test signal to modify the system parameters. 

( 3 ) Those that use some specific system interrelationship, such as 
control gain associated with neutral stability of one of the system 
components, to adjust the system parameters. 

Specific examples of each type can be found in reference (l). 

An automatic longitudinal-flight-control system that affords the 
pilot a constant ratio of pitch rate per pilct input co mma nd has been 
shown to be desirable (ref. 2). This ratio is easily obtained by a 
feedback control system with an integration in the forward loop. However, 
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it is necessary now to be concerned with the stability of the system, 
which can change due to changes in the controlled element. 

With the integration, the open-loop gain for a zero frequency exci- 
tation is infinite regardless of the variation in the controlled element; 
however, the system gain at a frequency other than zero is not infinite 
and is altered by the changes in the characteristics in the controlled 
system. One such characteristic is control effectiveness which varies 
with Mach number and dynamic pressure. This principle suggested the use 
of a system of the second type; that of measuring the vehicle response 
L to a high-frequency, small-amplitude, sinusoidal test signal producing 

1 oscillatory control-surface deflections. The control gain is automati- 

k cally varied to maintain the output pitch rate due to the test signal 

5 equal to a preselected constant. A narrow-band-pass filter is used to 

6 separate the response to the test signal from the response to other 
inputs. The purpose of this paper is to investigate the principle of 
operation; namely, to determine if adequate aircraft stability can be 
maintained by automatically varying the control gain level so that the 
controlled system response (as measured by vehicle pitch-rate response 
to a sinusoidal test signal) is kept constant. 

Included in this paper are the effects on system performance of 
variations in flight condition, system gain level, test-signal frequency, 
and sensor location. Results of an analog- computer study using a time- 
varying-coeff icient program for a flexible missile are presented to 
^ compare the control obtained with the adaptive system with that obtained 

with a fixed-gain system, and to determine the transient characteristics 
of the adaptive loop and the effect of these characteristics on the 
over-all stability of the flight-control system. 


SYMBOLS 


A(b) 


transfer function of airframe pitching velocity to control 
deflection 


a^ area of ith section, sq in. 

L 

^L.i lift coefficient of ith section, — — 

9 } 


dj_ normalized displacement of first bending mode at the force 

station of the ith section, in. /in. 
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error signal in adaptive loop 
filter transfer function 
force along Z-axis, lb 

hydraulic actuator transfer function 

displacement of reference point to which the first bending 
mode is normalized, in. 

adaptive -loop error shaping network transfer function 
moment of inertia about the Y-axis, lb-in. /sec 2 

imaginary portion of the complex variable s = a ± joo, 
radians/sec 

amplifier gain, dimensionless 

lift force of the ith section, positive upward, lb 

length between the missile center of gravity and the assumed 
force station of the ith section, in. 

pitching moment, about the Y-axis, in-lb 

mass, lb-sec^/in. 

2 / 

generalized mass of the first herding mode, lb-sec / in. 

forward-loop shaping network trarsfer function 
generalized force input to the first bending mode, lb 

dynamic pressure of air stream, Jb/sq in. 
dynamic pressure of exhaust stresm, lb/sq in. 

rectifier transfer function 
Laplace operator, per sec 

linear velocity component along the X-axis, in. /sec 
linear velocity along flight path, in. /sec 


CT\VJ1 4^ H 
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w linear velocity component along the Z-axis, in. /sec 

a angle of attack, radians 

local angle of attack of ith section measured at the assumed 
force station, radians 

5 control-surface deflection (tip control plus jet vane), radians 

€ forward-loop error signal 

£ damping ratio of linear second-order system particularized 

by the subscript 

9 angular displacement about the Y-axis, radians 

A ^ slope of normalized first bending mode at the force station 

of the ith section, per in, 

Aj slope of normalized first bending mode at station j, per in. 

a real portion of the complex variable s = a ± jo), sec"* 1 

t aircraft path time constant, sec 

0 angular displacement about the X-axis, radians 

i f angular displacement about the Z-axis, radians 

undamped natural frequency of second-order mode particularized 
by the subscript, radians/sec 

Subscripts: 

b first bending mode 

c input command 

G output measured at the center of gravity 

g output measured at sensor location 

r rectified version of signal 

s airframe short period 

t test signal 


px oenxng-vexucx by sxgnax ou one xnpux uesx-sxgnax ampxxxuae as a iuncxion 
of test- signal frequency are shown in figure 3 fo r the three flight 
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conditions of figure 2. This ratio Qg/& t will be referred to as 

closed- loop dynamic gain. The frequency and damping -ratio character- 
istics of the system components used in the calculation of figure 3 
represent those associated with the maximum and minimum control gains 
for that particular flight condition as determined by the root-locus 
study. Each value of closed- loop dynamic gain included between the 
boundary curves corresponds to a value of steady-state control gain 
associated with satisfactory short-period characteristics for that 
flight condition. 


Figure 4 is a superposition of the plots of figure 3 . The shaded L 
area starting at a frequency of about 8 radiai. s/sec represents combi- 1 
nations of test-signal frequency and closed-lcop dynamic gain common to 4 
aid the flight conditions. Thus, while there is no value of steady- 5 
state control gain common to the flight conditions investigated, there 6 


are certain test- signal-frequency ranges where there are values of 
closed-loop dynamic gain common to each, and each combination of fre- 
quency and gain represents a satisfactory steady-state control gain 
level. One such value of closed-loop dynamic gain at a particular fre- 
quency is used as a basis for the bias value and for the frequency of 
the test signal. Thus, as the control effectiveness varies, due to 
environmental change the control gain can be forced automatically to 
vary in such a manner as to maintain the closed-loop dynamic gain con- 
stant and in so doing will produce satisfactoiy short-period charac- 
teristics. 


Application to Missile Control 

The investigation was extended to include control of a flexible 
missile in order to further evaluate the system's ability to compensate 
for rapid changes in control effectiveness and to include the effects 
of flexibility and sensor location. A four-stage, solid-fuel launch 
vehicle with a longitudinal pitch-rate control system incorporating 
tip control flaps located on horizontal stabilizing fins and jet vanes 
immersed in the exhaust stream was chosen for the study. The resultant 
control effectiveness decreases rapidly after first-stage burnout due to 
the loss in effectivensss of the exhaust-stream vanes. A transfer func- 
tion relating missile pitching-velocity response to control-surface 
deflection was obtained by the method described in appendix A. Numerical 
values of the coefficients of the variables in the equations of motion 
were calculated as a function of trajectory flight time by a digital- 
computer program from unpublished vibration ani static-stability data 
for the test missile. 

The transfer-function values used in obtaining the system closed- 
loop equation are listed in table 3. The control system used in this 
portion of the study is still that shown in figure 1. 
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Root-locus plots showing the variations in the short -period and 
first-bending-mode characteristics as a function of control gain are 
shown in figure 5* Figures 5( a )> 5(t>), and. 5(c) represent three time 
periods along the trajectory for the same missile characteristics; 
whereas figure 5(d) is for a case with a higher frequency first bending 
mode, and 5(e) is for a case with a different sensor location. The 
sensor station numbers represent the number of inches between the sensor 
location and a reference point, which for the missile used, was the 
first-stage nozzle exit. Station 1+00 represents a location in the 
transition section between the first and second stages, and station 646 
represents a location in the section between the second and third stages 
The first-bending -mode slope value Aj at station 400 was 0.000 and 

was 0.0057 per inch at station 646. 

Since missiles generally operate with lower damping than aircraft, 
the limits of satisfactory short -period characteristics were extended 
from those associated with the aircraft to include damping ratios of 
0.2. The satisfactory region is shown in figure 5. A further condition 
imposed, was that the damping ratio of the bending mode should not be 
lowered below the missile’s basic value of 0.01 in attaining increased 
short-period damping, since the damping was initially low. 

The maximum and minimum values of steady-state control gain for 
the conditions shown in figure 5 are presented in table 4. No value 
of control gain is common to these flight conditions. It is of note, 
however, that a 54-percent change in first-mode bending frequency did 
not significantly alter the satisfactory control gain range although 
the value of control gain associated with neutral stability was changed. 
Also, a change In station location which represented a substantial 
change in the value of the bending-mode slope did not change the range 
of satisfactory values of control gain from those of the initial sensor 
location. 

This indicates that a margin of error in the knowledge of bending 
frequency and mode shape may be tolerable with this type of control; 
however, the magnitude of this error margin was not investigated further 
than the conditions shown. 

The closed-loop frequency response plots of the ratio of the meas- 
ured missile pitching-velocity-signal amplitude to the test-signal ampli 
tude as a function of the test-signal frequency are shown in figure 6 
for three flight conditions with the same missile conf iguration. The 
frequency and damping ratio characteristics of the system components 
used in the calculations were determined in the same manner as in the 
analysis of the longitudinal control of the high-performance research 
aircraft. Figure 7 is a superposition of figures 6 (a), 6 (b), and 6 (c). 
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The shaded region represents combinations of test-signal frequencies 
and closed-loop dynamic gains common to the various flight conditions. 

Thus, for example, if the adaptive loop w:>uld maintain the ratio 
of the output-to-input amplitude of a 15 radiais/sec test signal equal 
to -5 db, the missile short-period characteristics would be kept satis- 
factory, and the first bending mode would have increased damping through- 
out the atmospheric flight envelope. This -5 db value of dynamic gain 
assumes unity gain of the filter and rectifier, and if these gains were 
not unity, the bias value would have to be altered to account for the 
additional amount of gain present. 


ANALOG COMPUTER SIMULATION 


Description of Adaptive Loop 

The previous analysis methods were based on the assumption of 
linear system characteristics; however, any actual control system will 
exhibit some nonlinear characteristics due to such things as velocity 
limiting of the control actuator and position ...imiting of the control 
surfaces. An analog-computer simulation (see appendix B) of the missile - 
plus-control system analytically studied was conducted to indicate a more 
realistic estimate of the control available with the adaptive system, 
and to determine the transient characteristics of the adaptive loop. 

The transfer characteristics of the forward-loop shaping network N(s) 
and adaptive-loop filter F(s) and rectifier R(s) were realized by 
the use of d. c. amplifiers with suitable input and feedback networks. 
Three adaptive-loop band-pass filters of different values of selectivity 
were tested. The output -to- input relationship of the filters in Laplace 
notation is as follows: 


®t /®g = 


Ks 


2 2 


1 + 2^/cjjs + l / o > s 


(1) 


Suitable values of resistance and capacitance were chosen so that func- 
tions equivalent to damping ratios of 0 . 1 , 0 . 00 , and 0 . 005 , with natural 
frequencies of 15 radians/sec, could be simulated. The gain K was 
unity. In all further references the numbers 0.1, and so forth, on the 
filters refer to their equivalent damping ratio. The filter character- 
istics were such that filter 0 . 05 , while being more selective, had a 
slower transient response time by a factor of approximately as com- 
pared with filter 0.1. Filter 0.005 had such a slow response that it 
was unsatisfactory. 
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Ideally, the adaptive-loop shaping network l(s) should be an 
integrator with a very high gain so that the amplifier gain K v will 
be quickly adjusted to a value that will null the error signal e; how- 
ever, it was found necessary to add a signal proportional to the error 
itself to provide damping for the amplifier gain Ky. The gains of the 

two signals, e and J' edt, were adjusted to produce as rapid a change 

in amplifier gain per change in error signal e as possible and still 
not have excessive overshoot of the amplifier gain. These values were 
maintained throughout the remainder of the program. The resulting 
network transfer function is given by the following equation: 

e '/e = 0.5(l/s + 0.6) (2) 

The trajectory characteristics used in this study are shown in 
figure 8 and represent the portion of the flight between launch and 
second-stage ignition. 


Results of Computer Study 

Missile-trajectory time histories showing the response to a sharp- 
edge gust disturbance are shown in figure 9. All cases represent con- 
trol by the adaptive system with a sensor located at station 6 b 6 ; how- 
ever, figures 9(d.)^ 9( e )^ and 9(f) represent runs having increased 
sensitivity in the gain changing loop. Proper operation of the control 
system should cause the vehicle pitch rate to return to a value of zero 
following a disturbance. 

The disturbance was introduced into the simulation as a 1 ° step 
change in angle of attack; however, the notation on the figure is that 
of the more familiar horizontal wind velocity in ft/sec. This change 
in angle is related to horizontal wind velocity by the relationship 

^w = tan_1 T 2 ( 3 ) 

V v 

where V v is the missile vertical velocity component represented by the 
velocity along the flight path since the missile path is nearly vertical, 
U n is the horizontal wind velocity component, and An w is the change 

in angle of attack due to the wind. The 1 ° change is equivalent to 
wind velocities somewhat larger than the missile might be expected to 
encounter. 

Proper operation of the adaptive loop should maintain the rectified 
test signal 0£ r equal to the bias; however, this signal decreases 
somewhat during the later portion of the trajectory indicating a lag is 
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present in the adaptive loop. This lag is apparently caused by a limi- 
tation in the maximum rate of change of error signal e, allowing a lag 
in the change of the amplifier gain K^. However, the gain is increased 

sufficiently in order to damp adequately any pitching oscillation pro- 
duced by the disturbance, as shown by the trace. 

In order to try to increase the maximum 2 ate of change of e, both 
the bias and the rectified test signals were increased by a factor of 3 
in the simulation; thus the same reference ratio was maintained but 
the amplitude of the error signal e for a given change in the recti- 
fied test signal was increased. As shown in figures 9{&) , 9( e )> and 
9(f), this increase in maximum rate resulted in a more rapid increase 
in control gain but also appeared to cause an overshoot in the control 
gain, and secondary oscillations are evident in pitching velocity and 
bending traces. Thus, this increase resulted in the same conditions as 
previously mentioned in the discussion of the mechanization of l(s). 

The control gain was limited to a value of 80 in order to restrict the 
maximum amplitude of the control-surface oscillations to plus or minus 5° 
during the latter portions of the trajectory. 

The effect of filter characteristics can best be illustrated by a 
comparison of figures 9(^) and 9(h). The harmonic content of the dis- 
turbance signal was sufficient to cause an apparent increase in the 
test-signal amplitude as shown by the increase in 0^. r after 30 seconds 

of flight time had elapsed. This apparent increase was more pronounced 
in the 0.1 filter due to its wider bandwidth. This apparent increase 
in excitation signal amplitude caused the adaptive loop to lower the 
control gain. This decreased value was less than the minimum value for 
satisfactory stabilization, and following the disturbance increased 
amounts of pitch-rate and behding oscillations may be noticed. This 
apparent increase in test-signal amplitude appeared to be a function of 
both the selectivity and time constant of the filter. Some trade-off 
in the two is necessary in obtaining the best filter characteristics. 

Of the two shown, the 0.05 filter held the rectified signal 0^ r 

more nearly constant in the presence of extraneous signals and was felt 
to be a satisfactory compromise. 

Missile time histories showing vehicle response to a sharp-edge 
gust disturbance when controlled by a fixed-gain system are shown in 
figure 10. The fixed-gain type of system represents the type presently 
incorporated in the missile. A control-gain ^alue of 5 was used since 
it represented the value presently used in the* test missile. The sensor 
was located at station 646. 

The maximum amplitude of the pitch-rate and bending oscillations 
caused by the disturbance are similar for control by either system, but 
the time required to damp out the pitch-rate oscillations varies both 


cr\vn 4 ^ h \r* 
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with type of system and system components. The damping times for the 
various systems and combinations of components are compared in table 5- 
An evident advantage of the adaptive system is its ability to damp the 
pitch-rate oscillations in the latter portions of the trajectory. This 
increased damping will decrease the attitude errors present at the 
second- stage ignition. 

Another feature of the system is an ability to prevent instability 
induced by control sensitivity. The test-signal frequency needs to be 
near the frequency associated with neutral stability of some mode. As 
an increased value of forward-loop gain increases the tendency of the 
mode to become unstable, the dynamic gain should increase sharply. This 
increase should, in turn, call for a lower control gain through the 
adaptive loop. The operation of this feature was examined with the aid 
of the analog mechanization. The forward-loop gain is comprised of the 
network gain, amplifier gain, and actuator gain. A circuit was incor- 
porated in the mechanization to alter the network gain, which previously 
had been unity, without changing the amplifier and actuator gains. The 
network gain was given a step increase calculated to produce a value of 
forward-loop gain associated with instability, such as might be encoun- 
tered due to some partial electronic failure or apparent decrease in 
control effectiveness due to test-signal drift which would call for 
increased control gain through the adaptive loop. The results of this 
increase are shown in figure 11. Immediately following the network 
gain increase, initiated after 30 seconds of flight time, the vehicle 
became unstable (flat tops of the peaks represent the physical limita- 
tions of the recorder, not the response quantity itself); however, the 
amplifier gain (represented on the figure as control gain) decreased 
to its minimum value, and control was restored. Later, the rectifier 

signal GjI was less than the bias and called for an increase in con- 
x-, r 

trol gain. Since the amplifier gain was being multiplied by a large 
network gain, the forward-loop gain was very sensitive to amplifier gain 
changes, and after 4l seconds had elapsed, divergence started again but 
was quickly checked by a decrease in the amplifier gain. Network gain 
increases were initiated at several other points during the trajectories 
and the flight of figure 11 represents the most serious case; that is, 
the one that took the longest time to restore control. A similar run 
made with a fixed-gain system resulted in complete loss of control. 

The effect of sensor location on the adaptive control system T s 
performance is shown in figure 12. However, all the locations investi- 
gated had positive bending -mode slopes. It is evident that the bending 
slope at the sensor location has an appreciable effect on the magnitude 
of the bending induced by the test signal. The larger the bending 
slope at the sensor location the greater the induced bending. The maxi- 
m-urn amplitude of the excursions caused by the gust, and the time required 
to damp out these excursions, are similar, although the bias used in 
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all cases was based on a response ratio with a sensor location at sta- 
tion 646 . This similarity agrees with the results of the root-locus 
study which showed that the range of satisfactory control-gain values 
was similar for the different locations although the values associated 
with neutral stability were quite varied. 


CONCLUDING REMARKS 


The results of the investigation indicated that it is possible to L 

maintain satisfactory longitudinal short-period characteristics over a 1 

wide range of flight conditions, ranging from landing approach at a Mach 4 

number of 0,5 to a Mach number of 5.8 at an altitude of 150,000 feet 5 

for a high-performance research aircraft and over the atmospheric por- 6 

tion of a missile trajectory, by varying the steady-state control gain 
so that the closed-loop dynamic gain is kept constant. The dynamic 
gain is defined as the ratio of the vehicle pitch-rate response to a 


sinusoidal test signal at a particular frequency. Maintaining stability 
was possible despite the fact that common values of steady-state control 
gain associated with satisfactory longitudinal short-period character- 
istics did not exist for the wide ranges of flight conditions investi- 
gated because there were test-signal-frequency ranges where there were 
values of closed-loop dynamic gain common to each condition, and each 
value of dynamic gain represented a satisfactory value of steady-state 
control gain. 

The adaptive loop entailed only the addition to a rate-command 
control system of a logic circuit comprised oF a filter-rectifier net- 
work and a comparator-integrator network and vas simple in concept; 
however, some prior knowledge of the characteristics of a particular 
vehicle would be needed to determine the test-signal-frequency range 
and bias values that can be used for that vehicle. 

Also, care must be taken in locating the feedback sensor since the 
greater the value of the bending slope at the sensor location, the 
greater the induced bending. 

The sensitivity and selectivity of the f Liter somewhat limited the 
effectiveness of the system, but a suitable compromise was obtained and 
the adaptive system proved itself capable of closer path control of a 
missile during periods of low control effectiveness than fixed-gain 
systems yet did not overcontrol during periods of high control effec- 
tiveness. 
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In addition, the system demonstrated an ability to prevent insta- 
bility induced by increased control sensitivity due to some system 
malfunction. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., April 20 , 1961. 
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APPENDIX A 

DEVELOPMENT OF MISSILE PERTURBATION EQUATIONS OF MOTION 


The axis system employed and positive direction of forces, moments, 
and displacements are shown in figure 13 . 

The small perturbation equations for the rigid body longitudinal 
case can be written as follows (ref. 4): 

F z - *(* - V 0 9 q ) (Al) 

«Y * ¥ S 

The elastic equation for the first mode is (ref. 4) 

Q b = ^(h + 2^0^ + o>h 2 h) (A3) 

The force along the Z-axis is approximately equal in magnitude 
to the summation of the lift produced by the various sections of the 
missile and opposite in direction within the region of small angles of 
attach. 

The lift produced by each section can be expressed by: 

L i = qa i( C L a ). a i ( Al 0 


The aerodynamic force inputs were assumed to act at four force 
stations along the missile body plus at the centers of pressure of the 
horizontal fins, horizontal tip control surfaces, and the jet vanes. 
Substituting equation (A4) into equation (Al) yields 




(A5) 


The total pitching moment is equal to the summation of the products 
of the individual applied forces and their lever arms; therefore, equa- 
tion (A2) becomes: 


L 


h ltnvo 
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5 

Vo = ^ l a i^ C L a j + 

The generalized force associated with the first mode is composed 
of the summation of the amounts of work done by the individual forces 
in producing their relative displacements d^. Therefore, equation (Aj) 
becomes: 


qB6 (%) 6 16 + ^( S )-, 1 ' 


(A6) 


L 


5 

mb(h + 2^0^! + co^h) = ^>~ qa i( C L a ) ^ 

1=1 ' a/i 


+ 


n( c L s ) 6 d 6 



(AT) 


The local angle of attack can he expressed in the 
(ref. 4) : 


following manner 



\h 


(AS) 


further 


w 



a 


(A9) 


By applying equations (AS) and (A9) to equations (A5), (A6), and 
(A7), the following equations may be obtained: 


a = - 


q. 

mV. 


1+ ^g ai ( c 4^ 


q. 

mV 


L M74 di 

1=1 


h + 


i=l 



(A10) 
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Usually the equations will be solved numerically rather than in 
symbolic form with the various parameters combined to form a numerical 
coefficient for each variable. In order to aid in outlining the pro- 
cedure for solving these equations the coefficients will be represented 
by (n varies from 1 to 15 ). Equations fA10), (All), and (A12) 

may now be written as 

d = K-jO. + + Kjh + K^h + K^8 (A13) 

0 G = Kga + K ? 9 g + Kgh + Kh + K 1Q b (Al4) 

h * K-j_i a + + + J ** 1 + ^ 15 ^ (A15) 
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Equations (A13), (Al4), and (A15) may be written in Laplace trans- 
form notation 

(s - K^a(s) - K 2 0 g (s) - ^K 5 s + KjJh(s) = K^s) (A1 6) 

(s - K 7 )6 g (s) - Kga(s) - (Kgs + K^h(s) = K^Bfs) (A17) 

(s 2 - K 13 s - K l4 jh(s) - K lia (s) - K 12 9 G (s) = K^s) (Al8) 

The desired relationship is the transform of missile pitching 
velocity as measured by a sensor per control-surface deflection 
0g(s)/s(s). The sensing element is located at some station j on the 

body; thus the pitching velocity at this station is (from ref. b) 

0g = 9 q + (A19) 

The 0 (s)/b(s) relationship may be obtained from the following 
S 

relationship. 

0 g (s)/5(s) + Ajh(s)/5(s) = 9 g (s)/6(s) (A20) 

Equations (Al6), (A17)> and (Al8) may easily be solved by the method 
of determinants to obtain the desired relationships of equation (A20). 
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APPENDIX B 


Detailed wiring diagrams of the analog mechanization of the aero- 
dynamic equations and control- system simulation are shown in figures 14, 
15, and 16. In these figures, r indicates resistance measured in 
ohms; c, capacitance measured in pf ; and T h , hydraulic actuator time 

constant, measured in seconds. Initial conditions are indicated by IC. 


J H4 invo 
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TABLE 1.- CONTROL- SYSTEM TRANSFER FUNCTIONS FOR A 
HIGH-PERFORMANCE RESEARCH AIRPLANE 


N(s) 


O.^s + 1 H (s) = 1 

s( 0.05s + 1) ( 0 . 04 s + 1) 


Flight 

condition 


Mach 

number 

Altitude, 

ft 

A(s) 

0.5 

20, 000 

1 . 875 ( 3 - 51 s + 1) 
0 . 549 s 2 + 0 . 45 s + 1 

5-5 

70,000 

0.252(3. 91s + 1) 
0.0605s 2 +■ 0.0296s + 1 

5-8 

150, 000 

0.006(106. 5s + 1) 
0.628s 2 + 0.0119s + 1 


L 

1 

4 

5 

6 
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TABLE 2 .- MAXIMUM AND MINIMUM VALUES OF 
CONTROL GAIN OF A HIGH-PERFORMANCE 
RESEARCH AIRCRAFT 


Flight 

condition 

Maximum 

gain 

Minimum 

gain 

Mach 

number 

Altitude, 

ft 

0.5 

20, 000 

1.5 

0.2 

5-5 

70,000 

0.6 

0.4 

5.8 

150, 000 

10.0 

4.0 




TABLE 3.- CONTROL- SYSTEM TRANSFER FUNCTIONS OF A FLEXIBLE MISSILE 


2k 
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TABLE 4 .- MAXIMUM AND MINIMUM VALUES OF CONTROL GAIN 
OF A FLEXIBLE MISSILE 


Flight 

condition 

Sensor 

station 

radians /sec 

Control gain 

Mach 

number 

Altitude, 

ft 

Maximum 

Minimum 

1 

5,000 

646 

22.4 

11.0 

4.0 

3-25 

35, 000 

646 

22.4 

7.0 

2.0 

3-6 

100, 000 

646 

22.4 

200.0 

60.0 

3-25 

35,000 

646 

50.0 

6.0 

2.0 

3-6 

100, 000 

400 

22.4 

200.0 

60.0 
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Undamped 







Adaptive loop 
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function of control gain. 
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Figure 2.- Concluded- 
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Figure 3 -- Closed-loop frequency responses of adaptive pitch- rate aircraft control system at 
maximum and minimum acceptable control loop gain values for three flight conditions. 
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(c) Mach number, 0.5; altitude, 20,000 ft 
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Figure 4.- Frequency region of conmon ®t/®g ratios for adaptive aircraft control system as 
determined by superposition of frequency response plots of figures 3(a), 3(b), and 3(c). 
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number, 3.25; altitude, 35*000 ft; rate sensor station 646; (o fi * 1.9 radians/sec; 
0.08; = 22.4 radians/sec; ^ = 0.01 ; t & = 0.935 sec; and q = 31-4 lb/sq in. 
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Figure 7*- Frequency region of common 9^:0 ratios for adaptive missile control system as 
determined by superposition of frequency response plots of figures 6(a), 6(b), and 6(c). 
















Figure 8.- Missile trajectory characteristics 
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(a) Band-pass filter 0.05; disturbance, 21 ft/sec horizontal- velocity 
sharp-edge gust initiated at time 15 seconds. 


Figure 9* - Adaptive pitch-rate missile-control-system performance; tra- 
jectory portion from 5 seconds after launch to second-stage ignition 
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(d) Band-pass filter 0.05 and increased error sensitivity; disturbance 
32 ft/sec horizontal -velocity sharp-edge gust encountered at time 
20 seconds. 


Figure 9-- Continued. 
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(e) Band-pass filter 0.05 and. increased error sensitivity; disturbance, 
57 ft/ sec horizontal- velocity sharp-edge gust encountered at time 
30 seconds. 


Figure 9» - Continued 
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(g) Band-pass filter 0.1; disturbance, 21 ft/sec horizontal-velocity 
sharp-edge gust initiated at time 15 seconds. 


Figure 9* - Continued 
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Figure 10.- Fixed gain (b/e') pitch-rate missile-control-system performance, trajectory portion 

from 5 seconds after launch to second-stage ignition. 
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(b) Control gain 5; disturbance, sharp-edge gust of 57 ft/sec horizontal velocity initiated 

at time JO seconds; sensor located at station 646 . 

Figure 10.- Continued. 
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(c) Control gain 5; disturbance, sharp-edge gust of 59*5 ft /sec horizontal velocity initiated 

at time 55 seconds; sensor located at station 646. 
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(a) Filter 0.05; sensor at station 400 Aj = C. 000; disturbance, wind 

gust increasing from 0 to 21 ft/ sec horizontal velocity in 1 sec- 
ond initiated at time 50- 

Figure 12.- Trajectory time histories shoving the effect of sensor loca- 
tion on adaptive pitch- rate missile- control- system performance. 
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(b) Filter 0.05; sensor at station 646 Aj = 0.0037; disturbance, wind 
gust increasing from 0 to 21 ft/sec horizontal velocity in 1 second 
initiated at time 30 seconds. 


Figure 12.- Continued. 
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Figure 15. - Body-axes system employed with positive direction of forces, 

moments, and displacements shown. 







Figure 16.- Analog mechanization of adaptive -loop components 





